ANNALES 
UNIVERSITATIS MARIAE CURIE-SKŁODOWSKA 


LUBLIN — POLONIA 
VOL. XXVI/XXVII, 33 SECTIO AA 1971/1972 


Instytut Chemii UMCS 
Zakład Chemii Fizycznej 
Kierownik: prof. dr Andrzej Waksmundzki 


Władysław RUDZIŃSKI 


Jedno- i dwucząstkowe funkcje dystrybucji w ciekłych fazach 
adsorpcyjnych * 


OĄHO- u ĄByHACTMUHbIE (QYHKĄMM pPAaCHpeĄeJeHMA B XKMĄKUX AĄCOPÓNMOHHBIX 
cba3ax 


One- and Two-body Distribution Functions in Liquid Adsorbate Phases 


WSTĘP 


Ogólnie rzecz biorąc, opis teoretyczny cieczy różni się znacznie od 
opisu teoretycznego gazów z następujących powodów. Przede wszystkim 
brak tu tak wygodnego i powszechnie wykorzystywanego w teorii gazów 
małego parametru, jakim w przypadku gazów jest gęstość, nie można też 
zaniedbywać stanów, gdy większa liczba cząstek znajduje się w bliskiej 
od siebie odległości. Opis zaadsorbowanej fazy powierzchniowej, trakto- 
wanej jako ciecz, wydaje się więc bardziej skomplikowany niż opis adsorp- 
cyjnej fazy gazowej, choć jest on równocześnie nie mniej ważny. Należy 
bowiem zdać sobie sprawę z faktu, że wszystkie zjawiska adsorpcji wy- 
korzystywane w procesach wielkoprzemysłowych, to zjawiska adsorpcji 
przebiegające przy bardzo wysokich ciśnieniach wolnej fazy adsorptywu, 
a więc w takich warunkach, gdy faza powierzchniowa jest już cieczą. 

Dotychczasowy opis teoretyczny zjawisk adsorpcji koncentrował się 
głównie na tradycyjnym problemie wyznaczania izoterm adsorpcji, biorąc 
w miarę upływu lat do rozważań coraz bardziej skomplikowane układy 
w sensie wyrażenia na energię adsorbujących się molekuł. 

Był to więc opis globalny fazy adsorpcyjnej, bez wnikania szczegóło- 
wiej w sytuację wewnętrzną adsorbujących się molekuł. Jeśli już poja- 


* Publikacja niniejsza oparta jest na wynikach pracy doktorskiej autora, wy- 
konanej w Katedrze Chemii Teoretycznej UJ w Krakowie w r. 1970. 
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wiały się pewne stwierdzenia odnośnie tej sytuacji, to były to raczej arbi- 
tralne założenia modelowe niż wynik bardziej ogólnych rozważań. 

Celem tej pracy jest przebadanie własności ciekłych faz adsorpcyjnych 
na powierzchniach homogenicznych w oparciu o potencjałowy model zja- 
wiska adsorpcji i przy użyciu metod właściwych mechanice statystycz- 
nej. Badanie to polegać będzie na wyznaczeniu jedno- i dwucząstkowych 
funkcji dystrybucji, zestaw bowiem tych funkcji daje bardziej wyczerpu- 
jące informacje na temat sytuacji w fazach powierzchniowych niż jakie- 
kolwiek inne wielkości fizyczne. Wiadomo bowiem, że w takim przypad- 
ku, gdy energię oddziaływania międzycząsteczkowego można uważać za 
addytywnie złożoną z oddziaływań dwucząsteczkowych, cała termodyna- 
mika układu daje się wyrazić poprzez te dwie funkcje. Sama zaś postać 
tych funkcji jest źródłem informacji o walorach wybitnie poglądowych. 

Gdy chodzi o dotychczasowe prace na ten temat, to za niewątpliwie 
pionierską trzeba uznać pracę S. Ono [1] z uniwersytetu w Kyushu 
(Japonia). Metoda teoretyczna, którą się Ono posługuje, wywodzi się 
z jego wcześniejszych pokrewnych tematycznie prac. Faza powierzchnio- 
wa opisana jest za pomocą teorii dziurowej, z jednoczesnym zastosowa- 
niem przybliżenia Bragga-Wiliamsa, gdy chodzi o ocenę ilości par mole- 
kuła — dziura [2, 3]. Oddziaływanie z powierzchnią ciała stałego zostało 
z założenia ograniczone do pierwszej warstwy tego quasi-krystalicznego 
modelu, bezpośrednio sąsiadującej z powierzchnią. Ono zapowiedział 
wprawdzie rozszerzenie teorii na przypadek oddziaływań długozasięgo- 
wych, ale sądząc z piśmiennictwa i wyboru prac nadesłanych przez Ono, 
zapowiedź ta nie została na razie spełniona. Mankamentem teorii Ono 
jest też nieciągłość jednocząsteczkowej funkcji dystrybucji, a w pewnych 
przypadkach nawet jej niejednoznaczność. 

Temat ten został podjęty na nowo przed czterema laty w serii prac 
F. M. Kuniego [4,5] z uniwersytetu w Leningradzie. Punktem wyj- 
ścia jest w nich nowy sposób rozwinięcia korelacyjnych funkcji Ursella — 
Mayera w szereg funkcyjny, którego współczynniki rozwinięcia zależą od 
jedno- i dwucząstkowych funkcji dystrybucji w wolnej fazie adsorptywu. 
Metoda Kuniego, choć odznaczająca się dużą elegancją, ma jednak jedno 
poważne ograniczenie. Pozwala bowiem wyznaczyć jedynie asymptotycz- 
ną postać tych funkcji dystrybucji dla dużych odległości od powierzchni, 
takich że odwrotność tej odległości jest małym parametrem. Metoda jest 
zatem nieużyteczna w przypadku bliskich odległości od powierzchni, gdzie 
zachodzą zjawiska najciekawsze. 

Podjęte badania bazują na tunelowej teorii cieczy Barkera [6—8], 
obecnie najbardziej obiecującej teorii stanu ciekłego, a nadto na wyni- 
kach teoretycznych Barkera dotyczących układów pseudoliniowych w polu 
zewnętrznym. 
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Obecnie zostaną przedstawione zasadnicze cechy modelu tunelowego 
wraz z pryncypiami teorii Barkera dla układów liniowych w polu ze- 
wnętrznym, jako że stanowią one podstawy proponowanego w tej pracy 
sposobu. 

W modelu tunelowym układ dzieli się na pryzmy heksagonalne, z któ- 
rych każda zawiera jednakową ilość cząstek. Cząstki określonej pryzmy 
(tunelu) poruszają się niezależnie od położeń cząstek w pozostałych pryz- 
mach, ale tylko w obrębie własnej pryzmy. Ruch cząstek przypomina 
przeto ruch w tunelu. Zaniedbuje się też w modelu Barkera wymianę czą- 
stek między różnymi tunelami. Już te założenia pozwalają rozbić całkę 
konfiguracyjną na iloczyn jednakowych całek dla każdego tunelu z osob- 
na. Dalej następują szczegółowsze założenia. Zakłada się nadto, że oddzia- 
ływanie między cząstkami określonego tunelu jest takie, jak gdyby cząstki 
poruszały się wyłącznie po osi tunelu. 

Rozpatrując oddziaływanie określonej cząstki z cząstkami innych tuneli 
zakłada się możliwość ruchu tej cząstki w obrębie całego tunelu. Oddzia- 
ływanie jest jednak takie, jak gdyby cząstki poszczególnych serii tuneli- 
-sąsiadów były równomiernie ,,rozsmarowane” na ściankach cylindrów 
przechodzących przez osie tuneli należących do tych serii. Ilość tuneli w 
poszczególnych seriach i promień cylindra dla określonej serii dają się 
oczywiście wyznaczyć z geometrii podziału układu na pryzmy heksa- 
gonalne. 

Energia potencjalna oddziaływań daje się więc zapisać jako suma 
dwóch członów. Jednego, zależnego od odległości wzajemnej cząstek 
na osi tunelu, i drugiego, zależnego od odległości cząstek od osi tunelu. 
Matematycznie wyraża to możliwość dalszego rozbicia całki konfiguracyj- 
nej dla jednego tunelu na iloczyn M--1 całek, gdzie M jest ilością cząstek 
w tunelu. Wśród nich jest M jednakowych całek pojedynczych dla części 
płaskiej ruchu cząstek (w płaszczyźnie prostopadłej do osi tunelu) oraz 
jedna M-krotna całka, odnosząca się do ruchów podłużnych. 

Gdy chodzi o ruch płaski, to jest to, jak widać, ruch w polu o symetrii 
centralnej. Ów centralno-symetryczny potencjał wyznacza się sumując 
oddziaływania określonej cząstki w odległości R od osi tunelu, ze wszyst- 
kimi cylidrami określonymi dla tego tunelu. Poniżej zostaje przytoczony 
wynik Barkera dający wyrażenie v(R) na omawiany potencjał. 


V(R) = K jdz Jae U (|R + E — 2&,lcos¢ + 2? |)&! (1) 


> es 1 
gdzie: n;, l oraz §l to: numer serii sąsiadów, odległość między osiami tuneli 
oraz promień cylindra serii i. 

Z postaci potencjału widać, że omawiane całki jednokrotne, odnoszące 
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się do ruchu płaskiego, nie dadzą się wyrazić analitycznie, ale łatwo mogą 
być policzone w sposób maszynowy. M-krotna całka, odnosząca się do 
ruchów podłużnych, nie może być wyliczona w sposób bezpośredni na 
maszynach z uwagi na to, że krotność jej jest na ogół bardzo duża. Po- 
woduje to konieczność koncepcyjnego rozwiązania zagadnienia układów 
jednowymiarowych o dużej liczbie cząstek. Zagadnieniu temu poświęcone 
są dwie prace Barkera [9, 10], z których pierwsza dotyczy zwykłych 
układów pseudoliniowych bez obecności pola zewnętrznego, druga zaś — 
tych samych układów w polu. Przypadek pierwszy udało się Barkerowi 
sprowadzić do zagadnienia całki jednokrotnej. Z punktu widzenia autora 
niniejszej pracy bardziej interesujące są wyniki drugiej z omawianych 
prac Barkera. Właśnie przy okazji tej pracy rzucił Barker myśl o wyko- 
rzystaniu wyników w niej zawartych do zagadnień powierzchniowych. 
Jako pole zewnętrzne można bowiem w szczególności przyjąć pole sił 
adsorbentu. 

Należy jednak rozwiązać cały podstawowy problem, jak zagadnienie 
trójwymiarowej fazy powierzchniowej w polu zewnętrznym sprowadzić 
do punktu, w którym można by zastosować wyniki Barkera dla układów 
liniowych w polu zewnętrznym. 

Ściśle rzecz biorąc, układy rozpatrywane przez Barkera należy nazwać 
pseudoliniowymi, gdyż hamiltonian tych układów (część potencjalną), E, 
definiuje Barker następująco: 


a e — energia w polu zewnętrznym 
= 2e (r) + U(tyt + 1) U — energia oddziaływania dwu cząstek (2) 


gdzie, zgodnie z przyjętą w tej pracy notacją (brak strzałek), r; oznacza 
zespół trzech współrzędnych X;, y;, zi Rygor liniowości objawia się dopiero 
w zdefiniowaniu warunkowej całki konfiguracyjnej QN(r) w postaci: 


N e(r U(r, r p 

Qn (rss) =p] [exp [Oo 4 eha ar. dh; (3) 

Z1< Za. ZZNI 
Jest to całka konfiguracyjna układu, w którym jedna cząstka znajduje się 
stale w punkcie R 74,4, a pozostałe N cząstek przyjmuje położenia o skła- 
dowych zetowych w przedziale pomiędzy 0 i 2,44. Oczywiście taki zapis 
całki konfiguracyjnej wprowadza też założenie o istnieniu oddziaływań 
tylko między najbliższymi sąsiadami w tym pseudoliniowym układzie. 
Oznaczając dalej: 


FO) =exp| ER (a) 


H (r,r’) =exp pea (5) 


Jedno- i dwucząstkowe funkcje dystrybucji w ciekłych fazach adsorpcyjnych 271 


ze wzoru (3) otrzymuje się: 
Na» M2 N1(7)F()H (rr) dr (6) 


Definiując Q,(r) jako równe jedności, warunkową grandkanoniczną funk- 
cję rozdziału G(r) można zapisać w postaci: 


G= F NRO) (7) 
N=O 


Biorąc zaś pod uwagę równanie (6), otrzymuje sie następujące równanie 
całkowe dla tej warunkowej grandkanonicznej funkcji rozdziału: 
G()=1+2j)[G(r)F(r)H(r,r'”) dr” (8) 
ZZ ZZ 
Jest to warunkowa grandkanoniczna funkcja rozdziału dla układu rozcią- 
gajacego się w przedziale (0,z), gdzie z jest składową zetową cząstki 
o ustalonym położeniu. Zwykła grandkanoniczna funkcja rozdziału G* dla 
dyskutowanego układu, rozciągającego się w przedziale (o, Ż), będzie na 
podstawie równań (6—8) dana równaniem: 
G*=A ji F (r) G (r) Pr (9) 
ocz<Ż 
W celu wyznaczenia jedno- i dwucząstkowych funkcji dystrybucji defi- 
niuje się jeszcze warunkową grandkanoniczną funkcję rozdziału G'(r) dla 
układu rozciągającego się w obszarze (z, Ż). W analogii do równania (8) 
będzie mieć ona postać: 
G' (x,y,Ż—2)=1+A())G (x,y,Ż—z)F(x,y', Ż—z)H (r, 1’) d*r' (10) 
z<z'<Ż 
Wkład zatem stanów układu, w których jedna cząstka znajduje się na 
pewno w elemencie dr, otaczającym punkt r, do zwykłej grandkanonicz- 
nej funkcji rozdziału G* będzie: 
LF(r)G(r)G (r)d*r (11) 


Dzieląc równanie (11) przez G*, otrzymuje się prawdopodobieństwo zna- 
lezienia cząstki w elemencie dìr w otoczeniu punktu r. Stąd wynika wzór 
na jednocząstkową funkcję dystrybucji e(r). 

pO(r)=(G*)-1: A: F(r)G(r)G (r) (12) 


Podobne rozumowanie prowadzi do wzoru na dwucząstkową funkcję dy- 
strybucji oznaczaną dalej jako o*(14, Te). 


p” (rir) = (G*) "XG (r1) G' (re) F (r3) F (ra) H (r1 rs) (13) 


W dalszych rozważaniach przyjęto założenie, że model tunelowy Bar- 
kera jest właściwie odmianą teorii cel, o tyle doskonalszą, że dopuszcza 
ona możliwość znalezienia się bardzo wielu cząstek celi w jednym tylko 
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kierunku przy opisie ich wzajemnych oddziaływań. Stąd i z warunku mi- 
nimalizacji oddziaływań cząstek celi z resztą układu wywodzi się kształt 
tej celi jako pryzmy heksagonalnej. W zasadzie można by wyobrazić so- 
bie celę, której obszar stykałby się z powierzchnią adsorbentu, a ruch 
cząstek adsorptywu należących do tej celi odbywałby się po prostej pro- 
stopadłej do powierzchni przy opisie ich wzajemnych oddziaływań. W tym 
jednak przypadku nie można by absolutnie założyć jednakowego ,,roz- 
smarowania” cząstek na ścianach celi ani też jednakowej odległości ścian 
celi od osi cząstek. Wydaje się jednak bezsporne, zważywszy na kształt 
potencjału adsorpcji, jako funkcji odległości od powierzchni, że cela ta 
powinna mieć symetrię obrotową względem osi cząstek. Dalej odległość 
ścian celi od osi cząstek w pewnej od powierzchni odległości powinna być 
tym większa, im mniejsza jest gęstość adsorptywu na tej odległości od po- 
wierzchni. Na znacznych odległościach od powierzchni cela ta powinna 
przechodzić w normalny tunel, jak u Barkera. 

Nasuwa się od razu spostrzeżenie, że z uwagi na nieregularny kształt 
takiej celi, nie jest możliwy podział fazy powierzchniowej pomiędzy takie 
cele. Trudność ta okaże się jednak nieistotna. Dla znalezienia jedno- 
i dwucząstkowych funkcji dystrybucji wystarczy rozważać własności po- 
jedynczej celi. Biorąc powyższe argumenty pod uwagę, przyjęto, że cela 
ma kształt nieregularnego walca o promieniu równym W(z), gdzie W(z) 
jest wielomianem zmiennej z postaci: 


W(z)=l--az (14) 
l jest odległością pomiędzy sąsiednimi tunelami w cieczy Barkera o gęsto- 
ści Qo równej lim | POM) . Energia potencjalna cząstki znajdującej się na 


odległości z od powierzchni będzie sumą oddziaływań z powierzchnią i ze 
ścianami walca. Oznaczając przez e(z) potencjał adsorpcji, całkowitą ener- 
gię potencjalną cząstki w tunelu (w sensie pola zewnętrznego), oznaczaną 
dalej jako Vs, można zapisać w postaci: 


Vs(R, z)=v(R, z)+ e(z) (15) 
Powyżej V(R, z) jest potencjałem wytworzonym przez ściany celi. 


Ż 2% ut iziwsvUWYTIE AE WUW PN W 
vR =fah jagu |V Rte- tW- 2R e-h W cos? 1) x 


x p (2) W (16) 


Te założenia wskazują, że istnieje możliwość zastosowania wyników teo- 
retycznych Barkera. Cząstki tak zdefiniowanej celi stanowią bowiem 
układ pseudoliniowy w efektywnym polu zewnętrznym. Analizując od 
strony rachunkowej możliwości rozwiązania zagadnienia Barkera dla 
układów pseudoliniowych w polu zewnętrznym, trzeba stwierdzić, że 
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jest to zadanie dość trudne. Sprowadza się ono do rozwiązania równań 
całkowych i to pod warunkiem, że znane jest a priorż A. Znajomość 
A a priori można by założyć tylko w tym przypadku, jeśliby faza po- 
wierzchniowa miała duże wymiary w kierunku prostopadłym do po- 
wierzchni, wówczas bowiem można by utożsamiać A z jej wartością dla 
wolnej fazy adsorptywu. Dla dowolnego kształtu potencjału możliwe są 
tylko rozwiązania numeryczne, które stosunkowo prosto wyglądają jedy- 
nie dla równań z całkami jednokrotnymi. Jeśli już nawet samo rozwiązy- 
wanie równań całkowych nie jest łatwe, to trzeba dodać, że w propono- 
wanym modelu występuje szereg dodatkowych trudności. 

1. Biorąc pod uwagę fakt, że w wyrażeniu (16) na efektywne pole 
zewnętrzne V;3(R, z) występuje jednocząstkowa funkcja dystrybucji, ozna- 
cza to, że trzeba właściwie rozwiązać układ trzech równań całkowych (8), 
(10) i (12). 

2. W tym układzie równań całkowych prócz nieznanej wielkości A wy- 
stepuje jeszcze nieznana funkcja W(z). 

Biorąc to wszystko pod uwagę zaproponowano iteracyjną metodę roz- 
wiązywania tego układu równań całkowych, która analogicznie do metod 
chemii kwantowej może być określona jako metoda ,,samouzgodnionego 
pola molekularnego”. W zerowym przybliżeniu uważa się, że faza po- 
wierzchniowa zachowuje się jak gaz idealny, przeto zerowe przybliżenie 
p® (r) będzie mieć postać: 


p (r) =P, exp a (17) 


gdzie ə jest gęstością wolnej fazy adsorptywu nad powierzchnią adsor- 
bentu. Za zerowe przybliżenie A będzie się brać aktywność gazu idealnego 
o gęstości równej średniej gęstości fazy powierzchniowej traktowanej jak 
gaz idealny. 


Zz Wo (z) z Wo (z) Fo 
A) = |f az ar-e'"(@)-R| |f az | ar-r] (18) 


W,(z) jest zerowym przybliżeniem funkcji W(z). Znajdzie się go z warun- 
ku minimalizacji energii swobodnej A, danej oczywistym związkiem: 


A = 
pr NInA—1nG* (19) 


Energia swobodna jest funkcją a, ponieważ N oraz G* są funkcjami zmien- 
nej a. W zerowym przybliżeniu można położyć: 


pay er Z Wo(z) a) 
InG*=N,N,=)dzjdR:p_ (2)-R+27 (20) 
o o 


Minimalizując A ze względu na a znajduje sie W,(z). Nietrudno zauważyć, 


274 Władysław Rudziński 


że zerowe przybliżenie dla G(r) jest po prostu równe zerowemu przybli- 
żeniu dla G* dla układu zawartego w przedziale (o, z). 


2 wo (h) 
InG,(z)=)dh | dR-p® (h) R27 (21) 
o o 


Podobnie zerowe przybliżenie dla G'(r) jest równe zerowemu przybliże- 
niu G* dla części układu zawartej w przedziale (z, Ż), czyli: 


3 2 wo(h) 
InG,(2)=|dh | dRp$ (h)-R-2m (22) 
Wreszcie: 
F . (2) =exp| ra (23) 


Pierwsze przybliżenie er) znajdujemy następująco: 

1) przy pomocy F(z) znajduje się F;(r), używając g(z) oraz W,(z) do 
wyznaczenia V(R, 2); 

2) biorąc F;(r), G.(r),Go(r) w całkach wyrażeń (8) i (9) oraz h,, oblicza 
się G,(r), G’y(r) oraz Gy*; 

3) stosując (12) oraz G,(r), G'(r), Gy* znajduje się o (r); 

4) pierwsze przybliżenie 4 wyznacza się również z równania (12). 

Niech z, będzie taką odległością od powierzchni, na której faza po- 
wierzchniowa praktycznie ma takie same własności jak faza objętościowa 
cieczy Barkera. Najogólniej rzecz biorąc, mogą to być jednak odległości 
większe niż grubość rozważanej fazy powierzchniowej. W przypadku, gdy 
geometria porów adsorbentu dopuszcza możliwość, iż faza powierzchniowa 
może rozciągać się na odległość z, od powierzchni, sprawa jest stosunkowo 
prosta. Można bowiem założyć, że: 


p, = [W (Zo) “t ’ dR-P® (R, Zo)-R (24) 


a równanie (12) przekształcić do postaci: 


Pair |) Gy (r) Frl) r= 


V(Z,) 


(Zo) 
=w) | F,(R,z|1+2%, | Fy (r) G(r) H (r', R, zo)d'r']R-dR (25) 
o V(Zo) 


Wyżej przez V(z,) oznaczono obszar całkowania, stanowiący wnętrze celi 
aż do odległości z, od powierzchni adsorbentu. Oczywiście postać wyraże- 
nia (25) wynika z faktu, że przy warunkach z=Ż=z, wartość G'(z) równa 
się jedności. Po obliczeniu całek w (25) uzyskuje się równanie liniowe 
na À. Jeśli z powodów geometrii adsorbentu nie można założyć, że faza 
powierzchniowa dochodzi do odległości z, wówczas otrzymuje się nieco 
odmienny, choć podobny wynik: 
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p | GGMF,Mdr= 
V(Ż) 


W (Zo) 
=[W (z)] * | F,(R,Zy)|1+2, | G(r) Fr) H(r, R, z,)d,r |R:dr (26) 
S 1 15 (ż) 1 1 


Argumenty prowadzące do takiego rezultatu są następujące. Wiadomo, że 
jeśli nawet adsorpcyjna faza powierzchniowa może osiągać ograniczone 
tylko rozmiary w kierunku osi z, to musi istnieć pewien niewielki obszar, 
gdzie jej wymiary w kierunki osi z są praktycznie bardzo duże. Jest to 
obszar, przez który faza powierzchniowa kontaktuje się z wolną fazą 
adsorptywu, dzięki czemu zachodzi równowaga adsorpcyjna. Można tedy 
badać zachowanie się jednocząstkowej funkcji dystrybucji w tym właśnie 
obszarze i ustalić w zasadzie dowolny obszar zmienności tej funkcji. Na- 
leży jednak pamiętać, że znakomita większość oddziaływań międzycząst- 
kowych ma miejsce przy warunku ograniczoności układu. Dlatego całko- 
wanie należy przeprowadzić tylko po obszarze zajmowanym przez fazę 
powierzchniową. Tym samym funkcja G'(z) staje się i w tym przypadku 
równa jedności. 

Powracając do zagadnienia metody iteracyjnej, zarówno wyrażenie (25), 
jak i (26) stanowią liniowe równania na pierwsze przybliżenie A. Biorąc 
następnie o (r) i A, oraz dokonując tych samych kroków, jakie wymienio- 
no wyżej w punktach 1—4, itrzymuje się drugie przybliżenie na or) i A. 
Dalszy sposób postępowania jest już oczywisty. Na koniec należy je- 
szcze ostatnie wzięte przybliżenie na funkcję er) uśrednić po prze- 
kroju celi, gdyż z uwagi na kształt pola potencjalnego adsorpcji, funkcja 
ta powinna zależeć tylko od odległości od powierzchni. Warto wreszcie 
podkreślić, że w każdej następnej iteracji rachunki w oparciu o zespół 
grandkanoniczny odbywają się na bazie ciągle modyfikowanego układu 
w sensie przestrzennym. Jako ostateczny wynik, oprócz o%(z) oraz A, 
otrzymuje się też optymalny kształt celi. 


WYNIKI NUMERYCZNE I DYSKUSJA 


Gdy chodzi o wyniki numeryczne, ilustrujące dyskutowany wyżej 
model fazy powierzchniowej, to otrzymanie ich przy użyciu maszyny 
ODRA-1204 było możliwe dopiero po pewnym uproszczeniu modelu. 
Uproszczenie to polegało na: 

1) założeniu wyłącznie liniowej translacji cząstek celi przy obliczaniu 
potencjału oddziaływania ze ścianami celi; 

2) przyjęciu jednakowego przekroju celi w kierunku osi z. 

Gdy chodzi o warunek pierwszy, to wymuszony on został dość ogra- 
niczoną pamięcią maszyny. Chodzi tu głównie o zapamiętywanie tablic 
wartości funkcji H(r, r) z kolejnych iteracji. Warunek drugi wynikł z ko- 
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nieczności ograniczenia czasu liczenia do realnych wymiarów. Każdorazo- 
we optymalizowanie parametru a wydłużałoby czas pracy maszyny kil- 
kanaście do kilkudziesiąt razy. Przy takich uproszczeniach wyznaczenie 
e(z) trwa około 2 godz. Powyższa dyskusja kwestii technicznych ma nie 
tyle wyeksponować trudności związane z obliczeniami, ile fakt, że przy 
użyciu sprawniejszych maszyn obliczenia te mogą być prowadzone z wię- 
kszą i praktycznie zadowalającą dokładnością. Okazuje się przy tym, iż 
obliczenia numeryczne, nawet przy dyskutowanych wyżej ograniczeniach, 
prowadzą do ciekawych rezultatów, które obecnie zostaną przedstawione. 


U 


Ryc. 1. Wartości e(z) jako funkcja numeru iteracji 


Przede wszystkim należy zacząć od stwierdzenia, że proponowana pro- 
cedura iteracyjna jest zbieżna, przy czym różnice wyników między dzie- 
wiątą a dziesiątą iteracją są już tak niewielkie, że praktycznie można po- 
przestać na dziesiątej iteracji. Ilustruje to najlepiej ryc. 1. Należy dodać, 
że wykres dotyczy fazy powierzchniowej o grubości 5 średnic twardej 
kory oraz sytuacji, gdy gęstość wolnej fazy adsorptywu wynosi 0,05. Po- 


zostałe parametry: temperatura T=300°K, ĉo =1000°K dla potencjału 


oddziaływania z powierzchnią typu (3—9) oraz 2 


=400°K dla potencjatu 


oddzialywania miedzyczasteczkowego typu (6—12). 

Celem poznania najważniejszych cech jakościowych proponowanego 
modelu wykonano wiele obliczeń, między innymi by stwierdzić, jak zmie- 
nia się rozkład gęstości w fazie powierzchniowej w miarę wzrostu gęstości 
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wolnej fazy adsorptywu (ryc. 2). Wzajemne położenie krzywych na ryc. 2 
dowodzi prawidłowości przyjętego modelu. Zgodnie bowiem z faktem, że 
pole pod krzywą jest proporcjonalne do ilości zaadsorbowanych cząstek, 
uzyskuje się z obliczeń prawidłowy wynik, że adsorpcja rośnie w miarę 
wzrostu gęstości wolnej fazy adsorptywu. Stwierdzić ponadto można pew- 
ne uporządkowanie w fazie powierzchniowej, najbardziej wyraźne w bez- 
pośrednim sąsiedztwie powierzchni (pierwszy pik). 


? 


03 


Q2 


oO” 


1,0 4s 20 25 30 35 40 4S so z 


Ryc. 2. Przebieg funkcji o()(z) dla dwu gęstości wolnej fazy adsorptywu, w przy- 
padku gdy faza powierzchniowa osiąga grubość 5 średnic twardej kory 


Maksima funkcji ọ®(z) przypadają dla tych wartości z, dla których 
według szeroko rozpowszechnionych wyobrażeń o adsorpcji przypada po- 
łożenie „warstw adsorpcyjnych”, pierwszej, drugiej i trzeciej. Studiując 
ryc. 2 nietrudno zauważyć dość dziwne zachowanie się funkcji oz) 
w obszarze, gdzie z założenia kończy się faza powierzchniowa. Krzywa 
przybiera tam wartości poniżej gęstości wolnej fazy adsorptywu 0». Po- 
wód tego jest następujący. Biorąc fazę powierzchniową o grubości 5 śred- 
nie twardej kory, dokonuje się też maszynowego całkowania w tych gra- 
nicach, co jest równoznaczne z obcięciem oddziaływań w tym punkcie, 
gdzie z założenia kończy się faza powierzchniowa. Częściowo efekt tego 
obcięcia łagodzi ,,dowiazywanie” funkcji do wartości go na granicy fazy 
powierzchniowej w trakcie obliczania kolejnych przybliżeń aktywności A. 
Gdy z geometrii por adsorbentu wynika, że faza powierzchniowa może 
osiągać bardzo duże rozmiary w kierunku osi z, wówczas trzeba się zado- 
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wolić w miarę dokładnym wynikiem, gdyż z uwagi na czas obliczeń nie 
można dowolnie pogrubiać fazy powierzchniowej wziętej do rozważań. 
Z uwagi jednak na stosunkowo niewielki zasięg sił powierzchniowych na- 
leży się spodziewać, że już stosunkowo niewielkie grubości dadzą z prak- 
tycznego punktu widzenia zadowalające rezultaty. Istotnie, z rachunków 
przeprowadzonych dla gęstości @,=0,05 i trzech grubości fazy powierzch- 
niowej: 5, 6 i 7 średnic twardej kory, wysokości w maksimum pierwszego 
piku wynoszą: 0,26, 0,52 oraz 0,76. Widać stąd, że różnice wykazują ten- 
dencję spadkową. 

Analizując sytuację, gdy faza powierzchniowa może osiągać ograni- 
czone tylko rozmiary, przebadano, jak zmienia się rozkład gęstości wraz 
ze zmianą grubości fazy powierzchniowej. Zachowując przeto identyczne 
pozostałe warunki, rozszerzono obszar całkowania do 6 i 7 średnic twar- 
dej kory. Przebieg krzywych na ryc. 3 świadczy o znacznym wzroście 


10 20 30 40 sp 60 z 


Ryc. 3. Przebieg funkcji oW0)(z) gdy: 1 — faza powierzchniowa osiąga grubość 5 śred- 
nic twardej kory, 2 — faza powierzchniowa osiaga grubość 6 średnic twardej kory 


czystej gibbsowskiej adsorpcji przy wzroście grubości fazy powierzchnio- 
wej. Można też zaobserwować jeszcze odsunięcie się maksimów drugiego 
i trzeciego piku od powierzchni wraz ze wzrostem grubości fazy po- 
wierzchniowej. Jest to zapewne wynik zwiększonego oddziaływania mię- 
dzycząstkowego. Ponadto zauważyć można pewien wzrost uporządkowa- 
nia w przypadku grubszej fazy powierzchniowej. Wydaje się, że jest to 
rezultat pośredni wzrostu gęstości w fazie powierzchniowej. Jak z podanej 
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teorii wynika, w trakcie iteracyjnego wyznaczania 0 (z) uzyskuje się też 
kolejne przybliżenia na A. Podobnie jak w przypadku iteracyjnego wy- 
znaczania, kolejne przybliżenia dążą w oscylacyjny sposób do granicy 
z zadowalającą szybkością (ryc. 4). Powracając do jedno- i dwucząstko- 
wych powierzchniowych funkcji dystrybucji, należy stwierdzić, że pro- 
blem ich teoretycznego wyznaczania nabiera szczególnej wagi w sytuacji, 
gdy technika pomiarowa nie jest jeszcze w stanie dać żadnych wskazówek 
na ten temat. 


1 2 3 4 s 6 7 8 9 10 


Ryc. 4. Zmiany wartości A w kolejnych iteracjach dla fazy powierzchniowej 
o grubości 5 średnic twardej kory 
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PJAUNAWNH 
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PE3IOME 


Jims BbIBeĄCHMA YpAaBHEHM OJĘHO- M ĄBYdACTUAHOH (QYHKIMM pacrpe- 
JECHA B ZKMIKUX AĄCOPÓHMOHHBIX chasax B paboTe MCIHOJEZOBAHBI TEO- 
peTmueckue pe3yJIbTaTbI Bapkepa, OTHOCALIMECA K JIMHEHBIM CHCTE- 
MAM, a TAKKE K IIJIOTHBIM JKMNKOCTAM. TeopeTMUECKME Ą|OCTUZKEHKA AB- 
Topa MokaSaHbI Ha HMCJIOBBIX pacueTax OFHOUACTMUHOM cbyHKIWM pac- 
mpeyenenna. 


SUMMARY 


Barker’s theoretical results concerning the linear systems and dense 
fluids were applied in order to develop the equations for one- and 
two-body distribution functions in liquid adsorbate phases. Numerical 
calculations of the one body distribution function were performed to 
illustrate the author’s theoretical results. 
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